Feshbach resonance scattering under cylindrical harmonic confinement 
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A problem of collisions of atoms with two-channel zero-range interaction in an atomic waveg- 
uide is solved by using of a renormalization procedure. A matching of the solution to a solution 
of the related one- dimensional problem leads to relation between the one-dimensional and three- 
dimensional scattering parameters. The scattering amplitude and bound states for the confined 
system demonstrate differences from the related free and one- dimensional systems. 
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Introduction 



Quasi-one-dimensional atomic gases have been real- 
ized recently in elongated atomic traps (see Refs. 0,13); 
two-dimensional optical lattices (see Ref. Q), atomic 
waveguides (see Ref. y and references therein), and 
atomic integrated optics devices (see Ref. Isl and refer- 
ences therein). These systems attract recently increased 
attention due to their possible applications to atomic 
interferometry, quantum measurements, and quantum 
computations. Ultracold atoms under tight cylindrical 
confinement could reach the "single-mode" , or quasi- 
one-dimensional regime, where only the ground state of 
transverse motion is significantly populated at the ther- 
mal equilibrium. An analysis of two-body collisions in 
this regime in Ref. |^ demonstrates that the center-of- 
mass motion can be separated in a case of harmonic 
confinement and a zero-range interaction between free 
atoms leads to an effective one-dimensional interaction 
between confined atoms. The one-dimensional interac- 
tion strength demonstrates resonant properties as a func- 
tion of the ratio of the transverse width and the elastic 
scattering length. This confinement induced resonance 
has been analyzed in Refs. A related problem 

of two atoms under three-dimensional harmonic confine- 
ment has been considered in Refs. 0113, in El El. The 
works laiaiilllllllllllll analyzed a case of Bose 
atoms, describing their interactions by a Fermi pseudopo- 
tential. A case of confined fermions has been considered 
in Refs. 14, 15, , 16, 17]. 

The previous works deal with a case of interatomic 
interaction involving a single channel. However, a possi- 
bility of scattering properties control due to the effect 
of Feshbach resonance (see Ref. [H) brings attention 
to multichannel problems. A Feshbach resonance can 
appear if the energy of collision in an open channel is 
close to the bound state energy in a closed channel. A 
description of non-diagonal elements of a potential ma- 
trix in multichannel problems by a Fermi pseudopotcntial 
seems to be problematic due to its non-hermitian form. 
The present work, following Ref. E3> considers Hermi- 
tian d function interactions. Such interactions lead do 
divergences, requiring a renormalization of scattering pa- 
rameters. Probably, a multichannel zero-range potential 



method of Ref. could be used here as an alternative 
approach. 

The work is outlined as follows. A renormalization 
procedure for two-channel scattering under cylindrical 
harmonic confinement is presented in Sec. Section 
ini describes two models of two-channel one-dimensional 
systems: the atom-molecule one and the two-state one. 
The results are discussed in Sec. IIIII They include rela- 
tion of one-dimensional scattering parameters to three- 
dimensional ones, analyzes of scattering amplitudes and 
bound states. A system of units in which Planck's con- 
stant is ?i = 1 is used below. 



I. RENORMALIZATION 

Following Ref. E3j consider two-channel scattering 
of atoms with 6 function interactions under the exter- 
nal harmonic confinement described by the potential 
14onf (r)- Close-coupled equations for the wavefunction 
of the open channel ipa (r) and the amplitude for the sys- 
tem to be in the closed channel ■(/;,„ have the form. 



+ VaS (r) + Konf (r) 



m 

+ VamS (r) i; 

m 

E^Pm^DsDi^^m + K™^a (0) 



V'a (r) 



(1) 



Here m is the mass of the atom, E and r are, respectively, 
the energy and coordinate vector of the relative motion, 
Va is the strength of the open channel potential, Vam is 
the coupling strength, and D30 is a bound state energy 
in the closed channel. All the energies here are counted 
from the open channel threshold. Elimination of tpm from 
Eqs. Q leads to the one-channel Schroedinger equation 



E^Ja (r) 



-V^ + V,siE)6{r) + V,onfir) 
m 



■>Pa (r) 



(2) 

with the energy-dependent effective interaction strength 

\Vam\' 



Kff (E) = Va 



D 



3D 



(3) 



Let us at first reproduce, with some modifications, a 
renormalization procedure for collisions in free space 
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(Konf = 0) realized in Ref. [T3- The wavefunction in 
the momentum representation 



(q) = (27r) 
can be represented as 

V'a(q) = ^(q-po)+ 



-3/2 



d'^ripa (r) exp (— iqr) 



(4) 



{2n)-'Tfro.{po), (5) 



where po = \pmE is the colhsion momentum and the 
T matrix Tfrcc (po) obeys the Lippmann-Schwinger equa- 
tion 



1 + ^^fr- (^^o) / 



mq^dq 

pI- 



(6) 

The (5 function potential in the coordinate representa- 
tion leads to a constant potential in the momentum rep- 
resentation and, therefore, to a divergent integral in a 
Lippmann-Schwinger equation. This integral is regular- 
ized in Eq. by introduction of a momentum cutoff Pc- 
The zero-energy limit of Tftoc (po) should reproduce the 
dependence of the elastic scattering length on the exter- 
nal magnetic field B (see Ref. 



Jfroc (0) = — ( 1 + p ^ p 



(7) 



where a^u is the background elastic scattering length, 
-Bo is the resonant value of the magnetic field, and A 
is the phenomcnological resonance strength. The bound 
state energy D30 is proportional to the magnetic field, 
D3D — /i-B+const, where /i is the difference between the 
magnetic momenta of an atomic pair in the open and 
closed channels. As a result the non-renormalized pa- 
rameters in Eq. |^ can be related to aau, Bq, A, and /i 
as 



Va = — asD 1 asDPc 

TO V TT 



iKmP = — a3D/iA ( f - -asDPc 

TO \ n 



(8) 



D-sD = ^J' 



B-Bo-A + A[l - -a^DPc 

TT 



reproducing result of Ref. [l9j . In the limit pc ^ 00 the 
T matrix can be expressed as 



An 

Tfrcc (PO) = -« 

TO 

pI - mDsD 



asDPo - iPo - masnDsDPo + im {D3D + ^J.A) 



(9) 



Consider now a case of cylindrical harmonic confinement 
with the transverse frequency uj±, 



2 „2 



(10) 



where = x'^ + y'^ and x, y, and z are the components of 
a vector r. The Schroedinger equation can be written 
out as 



E^a (r) 

where 
Hi 



^-^ + H± + V,s{E)5{v) 



Id 1 92 



i'a (r) 



(11) 



TO \ op p op p ot) 



is the transverse Hamiltonian and 9 is the angular coor- 
dinate in the xy plane. 

The eigenstates of H±, denoted as IA^to^) (see Ref. 
i), satisfy 



\Nm,) = (iV + |m, I + 1) LU^\Nm,), 



(13) 



where to^ is the angular momentum. The cigcnfunctions 
can be represented in terms of generalized Lagguere poly- 
nomials. The value of jA^TO^) at the origin is zero for odd 
A^ and independent of A for even A^ = 2n, 



where 



\2nm.) = 



a± 



1 



/7ra_L 



(14) 



(15) 



is the transverse harmonic oscillator length. Therefore 
matrix elements of the interatomic interaction 



{2n'm'^\V,sS {r)\2n,m. 



1 



■VcsSomJom'Sz (16) 



are independent of n and n' . 

Let us represent the wavefunction -00 (r) as 

oc 

oo „ 

(r) = (27r)-i/' / (9) e''l2n0>. 

n 



(17) 



(The components proportional to [A^to^) with odd A^ or 
rriz 7^ are uncoupled and, therefore, could be excluded.) 
The coefficients V'n (q) satisfy the set of coupled equations 



— + {2n+l)u± 



00 

.-1-^V,S{E)J2 I dq'^n'iq') 



(18) 



For a collision of two atoms in a transverse state n with 
a relative axial momentum p„ = ^Jm[E ~ (2n + 1) uj±] 
the coefficients can be expressed as 



i'n' (q) ^ S{q-pn) Sn'n + 



TO 1 

g2 +i0 27r 



Tn'n (Pq) 



(19) 
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The elements of T matrix T„/„ (po) satisfy the Lippmann- 
Schwinger equations 



Tn'n (Po) 



1 



■V,s (E) 



i 1 

1 - y]! — (po) 

n"=0 ^ 



(20) 

where the level cutoff regularizes the divergent series. 
The right-hand-side of Eq. (|20|l is independent of n and 
n' . Therefore r„/„ (po) is independent of n and n' as well 
and has a form, 



Tn'n (Po) = Tconf (Po) 



n"— 



(21) 



A substitution of Eqs. (Q, ((HJ, and allows us to ex- 
press Tconf (po) in terms of the physical parameters as, 



TconfiPo) 



47ra I 



+ Y2[n- {poa±/2) 



ma± 1 mTfrco (po) 

1/2 2 



n=0 



-PcO-A 



(22) 



This expression has a finite limit at ric —> cxD,pc oo for 
a±_Pc — T^\/ Tic ~ (po^-l/S) ^- This limit can be written 
out as, 



Tconf (Po) 



where 



mTfree(po) V2' V 2 J 



(23) 



C ( I = lim 



g(n + a)"^/'-2(n, + a)i/' 



(24) 

with — 27r < argfn -I- a) < 0, is the Hurwitz zeta func- 
tion (see Refs. [Ml^)- Equation H23() is similar to Eq. 
(10) in Ref. 0, but the elastic scattering length is re- 
placed by an energy-dependent function ^Tfreo (po)- For 
a case of non-resonant scattering (A = 0) Eq. (|23|) re- 
produces the results of Refs. 0,I3- This derivation can 
also be implemented for a case of three-dimensional har- 
monic confinement, giving some justification to the use 
of an energy-dependent resonant scattering length in Ref. 
[T^ , although a Fermi pseudopotential is used there 
rather then 5 function. 



II. FESHBACH RESONANCE SCATTERING IN 
ONE DIMENSION 

A. Atom-molecule model 

Let us treat a bound state of an atomic pair in 
the closed channel as a molecule. A many-body one- 
dimensional system of coupled atoms and two-atom 



molecules can be described by the Hamiltonian similar 
to one used in a related three-dimensional problem (see 
Refs. El El), 



OO 

Ham= J dxl^¥^ (a 



Am dx^ 



Did {x) + 



(x) 

+ Ig^^l^ (x) *a (x) (x) + h.C 



(25) 



Here {y) and (a;) are the annihilation operators 
for the atomic and molecular fields, respectively, Ua is 
the strength of an interatomic interaction, and Did is 
the energy of the molecular state counted from the open 
channel threshold. 

A state vector of the two-atom system can be repre- 
sented as a superposition of atomic and molecular states, 

=e'^^[^r^L ix) 

OO 

+ i= J dy^o (y) {x - I) + I)] |o), (26) 

where P is a center-of-mass momentum, y is an inter- 
atomic distance, and |0) is the physical vacuum state. 

A substitution of Eq. 1)26(1 into the Schroedinger equa- 
tion 



i_ 

4m 



E I I*™) 



Harnmn 



(27) 



gives the coupled equations for the coefficients ipQ (y) and 



Evo (y) 



1 (Pipo 
m dy"^ 



Uaifio (0) + V25 vH -5 (y) 



(28) 



E^r = DiD^r + ^9Vo (0) 



Elimination of (p"™ from these equations leads to a single 
Schroedinger equation 



Evo (y) 



1 d'^(po 



C/cff (E) S (y) cpo (0) (29) 



m dy^ 

with the energy-dependent effective interaction strength 

2|.9P 



U,s{E) = Ua + 



E-D 



ID 



(30) 



B. Two-state model 



An alternative description of Feshbach resonances in 
one-dimensional systems involves atoms with two inter- 
nal states, a and 6, associated to the annihilation oper- 
ators ^'a (x) and ^'b (x) , respectively. The many-body 
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Hamiltonian of this system has the form, 

oo 



+Db¥, (x) (x) 



ct—a.b 
Ua 



*I (X) *I [X] (X) {X) 



+ Ubi'i (x) i'l (x) (x) ^-6 (X) 

{x)+h.C. 



(31) 



where Db is the energy mismatch between the states a 
and b and Uab describes transitions between the states 
on atomic cohisions. A state vector of two-atom system 
can be represented as a superposition of different atomic 
states 



I* 



ab\ 



dy 



|0), (32) 



where the coefficients (po (y) and ipf' (y) satisfy the set of 
coupled equations 



1 dVo 



TO dy^ 
+f/>f (0) 6iy) 

1 dVf 



Ua^O (0) 



(33) 



+ DiDipT (y) 



TO dy^ 

+ [Uab^o (0) + Ubpf (0)] S (y) (34) 



A Feshbach resonance can appear if one of the channels, 
ab, related to the coefficient ipf' (y), has a bound state. 
Such state can appear when Ub < and has the energy 
Did = Db - fU^. In following Db and \Ub\ wiU be 
tended to infinity keeping a fixed value of -Did. This 
limit allows us to neglect the third channel, bb, which 
will be infinitely distanced from the other ones. 

The coefficient (fl'^ (y) can be expressed from Eq. H34() 

as 

Uab 



2Kb/l 



-ipo (0) exp {-Kby) : 



(35) 



where Kb = \J m {Db — E). 

A substitution of Eq. 1)35(1 into Eq. H33|l leads again 
to the one-channel equation (|29|l . but now the effective 
strength is expressed as 



U,siE) = Ua- 



\Uab\^ 



2Kb/m + Ub 



This expression tends to Eq. (|30|l in the limit Ub 
while 



Db = DiD + jUl 



Uab^ 



VMUbi 



9- 



(36) 



-oo. 



(37) 



C. T matrix 



Equation (|29() has a simple exact solution. The wave- 
function in the momentum representation can be ex- 
pressed as 



ipo (q) = (27r) 



-1/2 



3 (pq -<i) + 



dyi^o (y) exp {~iqy) 
m 1 



- g2 + io 27r 



TiD (po) (38) 



where po — yrriE is the collision momentum and the one- 
dimensional T matrix satisfy the Lippmann-Schwinger 
equation. 



Tyd (po) = C/cff {E) 



1 + ^^1^ (Po) 



mdq 



.2 _ ,2 4_ jO 



Po - 9 + ' 



Its solution has the form 
TiD (po) = U,s (E) 



im , , 
2po 



(39) 
(40) 



III. RESULTS AND DISCUSSION 

A. One-dimensional parameters 

The limit of small collision momentum po corresponds 
to small values of the second argument of the Hurwitz 
zeta function in Eq. (|23|l . A substitution of expansion 

d\A -^^-C + 0{a),^/^\^-i^\, (41) 



where C ~ 1.4603 is the Olshanii constant (see Ref. '^), 
into Eq. H23|) leads to an expression which coincides with 
Eq. 1301) if 



U,s{E) 



na , 



-Tfrcc {E) 



1 



Cm 
A'Ka I 



Tfrco {E) 



(42) 



This expression is nothing else but the Olshanii formula 
(see Ref. 0) where the elastic scattering length and 
the one-dimensional interaction strength are replaced by 
^Tfrco (Po) and the energy-dependent strength Ucff{E), 
respectively. 

Equation H42|) . in a combination with Eqs. ^ and l|30|) . 
allows to relate the one-dimensional parameters Ua, Dio, 
and g to the physical parameters aso, A, Bq, and fi as 



Ua 



4a3_D 



l-C 



0-30 



2 2a3DfJ-^ 

Wl = — — 2— 



l-C 



(43) 



(44) 



Did ^^i{B-Bo)-uj^+ C^fiA 1 - (45) 
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The confinement-induced resonance (see Refs. HQ) not 
only scales here the interaction parameters, but also 
shifts the Feshbach resonance by the last term in Eq. 
gSJ. The term cj^ in Eq. gSJ reflects the shift of the 
continuum threshold due to transverse oscillations. 



B. Scattering 

Let us introduce dimensionless variables (the scatter- 
ing momentum k, the elastic scattering length a, the de- 
tuning 6, and the resonance strength d) as, 



Poa± 



E 



1 



.1 X 10 



.3 / EjnK) 1 
Y cj±(Mhz) 2 



b—ii 



B-Bo 1 



a=— w 2.8 X 10-^Jm(AMV)u±(MRz)a3D (nm) 



2aj_ti>j 



1.2 X 10 



_^a3rf(nm)^(^B) A(G)^m(AMU) 
v/wi(MHz) 



For example, for lo^ = 27r x 80 Khz, we have a k, 
2.3 X 10-3, h K. 64[B(G) - 853], and d k. 0.14 for 
the case of the 853 G resonance in Na; a « 0.07, 
& « 48 [B(G) - 1007], and d « 0.6 for the case of the 1007 
G resonance in ^^Rb; and a « 0.31, 6 « 39 [B(G) - 155], 
and d « 130 for the case of the resonance in ®^Rb. 

The T matrix for the scattering under cylindrical har- 
monic confinement (|23|l is expressed in terms of these 
dimensionless variables as 



afc — a5 -f d 



mw_L [afc2-a5 + d]C(l/2,-A:2) + fc2-6' 



(47) 



For small values of fc the zeta function can be approxi- 
mated by the expansion Eq. 1)4111 . leading to the expres- 
sion for one-dimensional T matrix [cf. Eq. (|40)l ] 



ifc3 - {ah + d)k 



(1 - Ca) fc3 + Mfc2 - [6 (1 - Ca) + Cd]k-i {ab - d) ' 

(48) 

A study of the opposite, high k, limit requires an expan- 
sion for ^(i,a;) at high values of a. For a > such 
expansion can be obtained from the integral representa- 
tion (1.10.7) in Ref. ^ 



1/2 



24 



-a 



-3/2 



(49) 



At a < the function C (5,0;) has complex values with 
the imaginary part represented by a finite sum (see also 



Ref. a). 



I - a 



;^(|a| -n) 



ri=0 



-1/2 



(50) 



where [\a\] denotes the integer part of \a\. In the limit 
|a| — > 00 this sum can be approximated by the Hurwith 
zeta function itself [see Eq. (|24|l ] . leading to 



(51) 



A substitution of the leading term of this expansion into 
Eq. llTr|) gives us 




ak"^ — ab + d 



muj± 2ak^ — ifc^ — 2 [ab — d) k + ib 



(52) 



This expression differs from Tfrcc (see Eq. ^) only by 
a factor due to different definitions of T matrix for the 
scattering of the confined and free atoms. 

The one-dimensional scattering can be characterized 
by reflection probability 



R — I /even I — I -^confl 

2po 



and by the phase of the scattering amplitude 



X = argTconf (fc) - -. 



(53) 



(54) 



A comparison of results obtained with the exact ex- 
pression (|47|l and the approximate ones Eq. H48|l and Eq. 
(|52|l is presented in Figs. mill and|21 As one can see, the 
one-dimensional approximation (|48|l is appropriate only 
at small values of k. The three-dimensional approxima- 
tion (|52l) reproduces an average behavior at high fc, but 
it cannot reproduce the jumps appearing due to opening 
of transverse channels at integer values of fc^. 



C. Bound states 

The bound state energies of two atoms in the atomic 
waveguide are given by the poles of Tconf (po) on the pos- 
itive imaginary axis of pq. Equation H47|) leads to the 
transcendent equation in x = — i/c 

-^^^^ = -c(l,A,x>0, (55) 
ax^ + ab — d \2 J 

determining the dimensionless binding energy 

e& = a;^, (56) 




FIG. 1: The reflection probability [(a), see Eq. I53ll 1 and the 
scattering phase [(b), see Eq. (1541 1 calculated as functions of 
the dimensionless collision momentum k [see Eq. I46ll 1 with 
the parameter values a — 0.1, 6 = 0, and d = 1. The solid, 
dashed, and dot-dashed lines are related, respectively, to the 
exact expression 1471 . to the one-dimensional approximation 
(1481 . and to the three-dimensional one 1521 . 




FIG. 2: The same as in Fig. Q but for the parameter values 
a = 0.1, b = 10, and d = 1. 




FIG. 3: The same as in Fig. Q but for the parameter values 
a = -0.1, b = 10, and d=l. 



such that the bound state energies are 



E : 



-2UJJ_ Eb 



1 



(57) 



For shallow bound states {\E ~ u!±\ ^ u!±, x ^ 1), a 
substitution of the expansion (|41|l gives a cubic equation 

(1 - Ca) x^ + ax^ + [b (1 - Ca) + Cd] x + ab~d = (58) 

for bound states in the related one-dimensional system. 
A similar equation has been used in Ref. ||2^ for evalua- 
tion of binding energies of "bosonic mesons" . 

For deep bound states {\E\ 3>w^, x ^ 1) a substitu- 
tion of the leading term of the expansion (|49|l into Eq. 
H55|l results in a cubic equation 



2ax-^ 



+ 2{ab-d)x-b^0, 



(59) 



determined bound states of free particles (similar equa- 
tion has been considered in Refs. 25] ). 

A general behavior of bound states can be seen from 
a qualitative analysis of Eq. (|55() . Its right hand side 
varies monotonically from — cx) to cx) for < x < 00. 
The left hand side is a monotonic function of x for x > 
unless b < d/a. Therefore Eq. (|55|) has two real posi- 
tive roots for b < d/a and only one such root otherwise. 
The two roots correspond to two bound states. One of 
them tends at 5 —> —00 to the bound state in the closed 
channel. At 5 ^ 00 it tends to the bound state of the 
related one-channel system analyzed in Ref. ■ The sec- 
ond bound state tends to the state of the one-channel 
system at b ^ 00 and vanishes at the continuum thresh- 
old at b = d/a. A numerical solution of Eq. (|55|l (see 
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-10 -5 5 

b 



FIG. 4: The binding energy [see Eq. I|56|l ] calculated as func- 
tions of the dimensionless detuning b [see Eq. 1461 1 with the 
parameter values a — 0.1 and d = 1. The solid, dashed, 
and dot-dashed lines are related, respectively, to the exact 
expression 1551 . to the one-dimensional approximation I58L 
and to the three-dimensional one 1591 . The liner plot (a) 
demonstrates a general behavior, while the logarithmic one 
(b) highlights the near-resonance region. The solid and dot- 
dashed lines almost coincide in the part (a). 

Figs. 01 and |SJ| confirms the results of the quantitative 
analysis, demonstrating an intermediate behavior of the 
confined problem between the one-dimensional and three- 
dimensional approximations. In a case of a positive back- 
ground scattering length (see Fig. 0)) the bound states 
of the confined problem tend at |6| — > cxd to the bound 
states of the three- dimensional model, while the one- 
dimensional one has a single bound state only. An op- 
posite situation takes place at negative scattering length 
(see Fig. ISJ, when the three-dimensional model has one 
bound state only, and the bound states of the confined 
problem tend at j&j — > oo to the two bound states of the 



one-dimensional model. 



Conclusions 

A problem of two-channel scattering of atoms under 
cylindrical harmonic confinement can be solved using a 
renormalization procedure. 

Many-body problems involving two-channel scatter- 
ing can be described in one dimension by two models: 
the atom-molecule one and the two-state one. Parame- 
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FIG. 5: The same as in Fig. 0] but for the parameter values 
a — —0.1 and d = 1. The lower solid and dashed lines almost 
coincide in the part (b). 

ters of these models can be expressed in terms of three- 
dimensional scattering parameters, using the solution of 
the confined problem. Scattering amplitudes and bound 
states of the confined system incorporate both proprieties 
of the related one-dimensional and three-dimensional sys- 
tems, as well as specific peculiarities. 
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